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1.  INTRODUCTION 


The  generation  of  neutral  excited  state  atoms  or  molecules  in  a  liquid  helium  bath  via  collisions 
with  alpha  particles  was  initially  reported  by  Surko  and  Rcif  (1968).  Subsequent  experiments  utilizing 
discharges  from  beta  emitters,  again,  submerged  in  liquid  helium,  also  found  a  neutral  entity  in  a  long- 
lived  excited  state  (Rayfield  1969;  Mitchell  and  Rayficld  1971).  This  excited  atom  or  molecule 
produced  a  Uet  ion  and  an  electron  at  the  liquid  surface.  It  was  suggested  that  tliis  excited  species 
was  either  the  helium  2^S  atomic  state  or  the  a^Z*  diatomic  state  which  is  known  to  be  bound. 

Calvani,  et  al.  (1972),  generating  the  neutral  entities  from  an  alpha  source,  set  a  lower  limit  of  0.1  sec 
on  its  natural  lifetime  (t). 

A  more  recent  experimental  study  by  Melirotra,  Mann,  and  Dahm  (1979),  concluded  that  the 
neutral  excited  species  was  the  a^Zu  molecular  state,  i.c.,  the  lowest  energy  excited  state  in  Hc2,  and 
not  the  2^S  atomic  state.  These  workers  predict  a  lower  bound  on  the  lifetime  of  10  sec  in  liquid 
helium.  Tliis  lifetime  supports  die  assignment  of  the  electronic  state  to  the  molecular  a^Zu  rather  than 
the  atomic  2^S  state  wliich  has  an  experimentally  known  lifetime  of  15  p.sec  in  the  liquid  helium 
(Mchrotra,  Mann,  and  Dahm  1979). 

Another  interesting  aspect  of  the  excited  molecule  is  that  is  reportedly  forms  a  microscopic  bubble 
in  the  liquid  (Dennis  ct  al.  1969;  Hickman  and  Lane  1971)  with  a  theoretically  estimated  diameter  of 
12.5  A  (Hansen  and  Pollock  1972).  The  size  of  the  bubble  is  attributed  to  the  size  of  the  2s  Rydberg 
orbital  comprising  the  He  2^S  stale  (Gubennan  and  Goddard  1975). 

In  view  of  the  large  difference  (a  factor  of  100)  predicted  by  the  two  experiments  (Calv;ini, 
Maraviglia,  and  Mc.ssana  1972;  Calvani  et  al.  1974;  Mchrotra,  Mann,  and  Dahm  1979)  for  the  lower 
limit  of  the  lifetime  of  the  a^Z^  slate  in  liquid  helium,  high  quality  ah  initio  calculations  were 
undertaken  in  an  effort  to  clarify  this  situation. 

In  this  study,  the  lifetime,  x,  for  the  spin-forbidden  transition  a^Z^  — >  X’Z^  is  obtained  from 
calculations  employing  state  averaged  mulliconfiguration  SCF  (MCSCF)  plus  configuration  interaction 
(Cl)  wavcfunctions  to  describe  the  appropriate  zeroth-order  states.  In  order  to  calculate  tliis  spin- 
forbidden  lifetime,  the  spin-orbit  (S-0)  induced  perturbation  T*'  to  each  zcrotli-order  state  will  be 
calculated  using  the  full  microscopic  Breil-Pauli  Hamiltonian  (Bethe  and  Salpelcr  1977).  A  recently 
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implemented  method  (Yarkony  1986,  1987)  which  employs  the  symbolic  matrix  element  approach  of 
Liu  and  Yoshiminc  (1981)  in  the  evaluation  of  the  S-0  matrix  elements  will  be  used  to  evaluate  the 
H*'  directly  from  a  system  of  linear  equations  in  the  configuration  state  function  (CSF)  basis.  This 
technique  has  been  used  successfully  for  studying  spin-forbidden  transitions  in  other  molecules 
(Yarkony  1986,  1987). 

2.  METHODS 

The  Cl  method  used  to  obtain  tlic  zeroth-order  wavefunctions  is  the  symbolic  matrix  element, 
dircct-CI  method  of  Liu  and  Yoshiminc  (1981).  The  molecular  orbitals  (MOs)  needed  as  a  basis  set 
for  the  Cl  expansions  were  obtained  from  a  .state-averaged  mulliconfiguration  self  consistent  field  (SA- 
MCSCF)  approach.  The  SA-MCSCF  procedure  is  the  general  second-order,  density  matrix  driven 
MCSCF  algorithm  of  Lengsfield  (1982).  From  the  SA-MCSCF  procedure  one  obtains  a  set  of 
molecular  orbitals  detemiined  by  minimizing  the  energy  functional 

where  the  4^^’s  are  the  eigenfunctions  of  11°,  the  non-relativistic  Hamiltonian  operator,  in  tlic  space  of 
the  MCSCF  expansion 

'v:  =  E  c.v. .  (2) 

I 


Above,  the  w^’s  are  the  non-negative  weighing  factors  for  the  electronic  states  which  do  not  vaiy'  as  a 
function  of  intemuclcar  .separation,  and  the  v,  arc  the  CSFs  composed  of  the  state-averaged,  optimized 
MOs.  The  weights  were  chosen  to  provide  a  balanced  de.scription  of  the  states  of  interest.  The 
sensitivity  of  the  multi-reference  Cl  rcsulLs  to  a  particular  choice  of  Wj^’s  in  the  MCSCF  optimization 
was  tested  by  varying  the  weights  and  by  comparing  our  results  to  both  experimentally  derived 
spectroscopic  parameters  and  to  spectroscopic  parameters  obtained  in  earlier  ab  initio  calculations 
which  employed  a  separate  MCSCF  procedure  for  each  state  (Konowalow  and  LcngsOcld  1987a, 
1987b). 

2.1  Spin-Orbit  Inicraciions.  The  spin-orbit  part,  H*®,  of  the  microscopic  Breit-Pauli  interaction  is 
given  by  (Bclhe  and  Salpeter  1977) 


= 


e 


(3) 


2mc^ 


■>/ 


f 

->>1 
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Pi 

‘j 
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V 
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Due  to  helium’s  small  nuclear  charge,  the  S-0  interactions  arc  expected  to  be  small  and,  therefore, 
well  suited  for  treatment  by  first  order  perturbation  theory'.  The  total  perturbed  wavcfunction  for  state 
I  is  given  by 

'F/  =  'p;  +  .  (4) 


with  4^  being  the  zcroih-order  wavcfunction.  The  usual  spectral  representation  for  the  first-order 
correction  due  to  S-0  effects  is 


M 


(e;  -  £/) 


4^; . 


(5) 


The  summation  over  the  L  electronic  slates  is,  in  principle,  infinite.  One  often  used  approach  to 
solving  for  4^1  is  to  calculate  explicitly  the  wavcfunciions  for  a  relatively  smtill  number  of  excited 
states  thereby  drastically  truncating  L,  This  might  cause  one  to  miss  important  contributions  to  4'', 
from  die  omitted  states. 


Within  a  given  CSF  space,  this  "omitted  states"  problem  is  eliminated  by  solving  for  4'j  directly 
from 

(W"  -  £)  4'',  =  -W*'’4^;  .  (6) 

Equation  6  can  be  Iransfonncd  into  matrix  foim  as 

(//"  -  E)V'  =-tr  C'  (7) 

where  it  must  be  emphasized  that  H’  and  }_f"  arc  matrices  with  elements  fonned  over  CSFs,  not  over 
eigenstates.  The  vectors  v'  and  C'  arc  defined  as  the  coefficients  for  the  first-  and  zeroth-order  parts 
of  T,: 

4^;  =  E  C!  y.(x-) 

=  E  Vj  V/K-)  . 


3 


(8u) 

(8b) 


The  K  and  k’  label  ihe  spatial  symmetries  to  which  the  CSFs  belong,  and  in  general,  k^k'.  Equation 
7  forms  a  large  set  of  linear  inhomogeneous  equations  which  arc  solved  to  obtain  V  by  a  variant  of  the 
method  suggested  by  Poplc,  ct  al  (1979). 

2.2  Perturbed  Wavofunctions.  I’lie  following  perturbations  to  ‘P°(X'Zgo+)  and  ‘P°(a^Zui)  arc 
calculated 


(9a) 


(9b) 


where  the  first-order  corrections  arise  from  the  S-0  interactions 


y‘(^n,^;A''E;„j  :  <  >  f2  =  o* 

4^‘('n /.^’E:)  :  <  a’E:/l//"l'n  ,  >  f2  =  1  . 


Where  tlic  quantum  number  Q  (Q  =  A  -t-  S^),  the  z-componcnl  of  the  total  orbital  and  spin  angular 
momentum  is  conserved.  Below,  the  first-order  wave-functions  w'ill  be  abbreviated  as  ^’('0^^^^  and 

4^'('n„). 


2.3  Electronic  Transition  Dipole  Moment.  In  order  to  calculate  the  lifetime  of  tlic  a^E^  — >  X'Eg 
transition,  the  electric  transition  dipole  moment  Pi(a^Eu.X'Ep.  defined  by 


p,(u'e:,a''e;)  =  <  'p  (cj^e:,)  (.y'e;j  > 


(10) 


is  required.  The  quantity,  is  the  shift  operator  fonn  of  the  total  electric  dipole  moment  operator 
which  has  components  Substituting  the  perturbation  expansion  for  each  state  in  Equation 

10  gives,  to  first-order. 
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(11) 


Since  the  lower  state  in  this  transition  is  largely  repulsive  (possessing  only  a  very  shallow  van  dcr 
Waals  well),  we  need  to  obtain  the  vibraiionally  averaged  transition  dipole  moment  between  a  bound 
electronic  state  (here  the  with  vibrational  wavcfunction  Xv>(l^)  ^>^^1  a  repulsive  state  (here  the 
X'Zg  state)  with  a  continuum  vibrational  wavcfunction  Xk'*(l^)  's 

5  V'  =  <  x..(/?)iMi(«'5::,A:‘E;)ix, .,(/?)  >  (^2) 


where  k”  represents  tlie  energy  for  the  continuum  state  (van  Dishocck,  Langhoff,  and  Dalgamo  1983; 
van  Dishocck  and  Dalgamo  1983).  Xk.  ”(^)  Xv’(l^)  ‘’re  obtained  by  numerically  solving  the  radial 

Schroedingcr  equation  for  nuclear  motion  while  ignoring  rotational  effects.  The  vibrational 
wavcfunction  for  the  bound  state  is  normalized  to  unity  and  the  continuum  wavcfunction  is  defined  by 


x*(/?)  = 


% 

[Kk^ 


s\n(kR-n) 


(13) 


where  n  is  a  phase  shift  factor  and  p  the  reduced  mass  of  He2. 

The  Einstein  coefficient  for  spontaneous  emission  from  the  v’  to  tlic  k”  vibrational  state  is  (van 
Dishocck,  Langhoff,  and  Dalgamo  1983;  van  Dishocck  and  Dalgamo  1983); 


=  (2.I419a:70"’)•A£’(m<)^5^,,/  . 


(14) 


The  radiative  lifetime  for  the  v'  level  is  obtained  by  integrating  Equation  14  over  k". 

3.  DETAILS  OF  CALCULATIONS 

The  Gaussian-type  basis  set  is  csscnlitdly  that  used  by  Sunil,  ct  al.  (1983)  in  an  earlier  theoretical 
study  on  the  excited  states  of  HC2  with  tw'o  exceptions.  A  single,  primitive  p  function  has  been  added 
with  its  exponent  optimized  in  increments  of  0.001  to  give  the  lowest  energy  for  the  F'fT^j  at  R  =  2.00 
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bohr  (i.c.,  near  Tj.).  The  Cl  pari  of  ihe  oplimizalion  used  MOs  obtained  from  a  MCSCF  calculation  on 
the  F'n^j  state  (i.c.,  no  state  averaging).  Tliis  additional  p  function  was  deemed  necessary  due  to  an 
unacceptably  large  AECF’n^  -  a^Xu)  at  R  =  2.00  bohr  when  compared  with  the  experimental 
between  these  two  states.  The  orbital  exponent  for  the  more  diffuse  d-function  was  also  changed  to  be 
consistent  with  a  basis  set  used  in  an  earlier  study  on  Hcj  conducted  in  this  laboratory  (Konowalow 
and  Lcngsficld  1987a,  1987b).  The  final  atomic  basis  set,  reproduced  in  Table  1,  consists  of 
(10s,6p,2d)  primitives  contracted  to  (7s,5p,2d),  for  a  total  of  34  atomic  basis  functions  per  atom. 

The  ctUculations  are  performed  in  symmetry  with  the  appropriate  averaging  of  states  in  the 
SA-MCSCF  to  give  wavefunctions  winch  transform  according  to  symmetry.  In  the  states 
transform  according  to  the  irreducible  representations  (IRREPs)  X'Sg('A^),  a^X^C^BiJ,  b^FIgCx; 
yi^Bjp,  and  F'O^Cxi’Bj^,  yi'B.J. 

The  SA-MCSCF  is  of  the  CAS  type  wherein  the  four  electrons  arc  distributed,  in  all  possible 
ways,  amongst  the  lowest  three  MOs  from  IRREPs  3^(0^,)  and  b,^(Oy),  and  the  lowest  MO  from 
*^2u(^uy)’  t-'^’nsistcnt  with  space  and  spin  symmetry'  restrictions.  The 

state  averaged  energy  is  then  optimized  according  to  Equation  1,  including  the  states  X’Xg.  a'X*. 
b'rig^,  b^n^y.  F'riy^,  and  F'n^y.  Two  different  weighing  schemes  were  used  in  this  study.  The 
weights  w=  (2,  2,  1,  I,  I,  1),  and  w  =  (1.5.  1.5,  1,  1,  1,  1)  were  employed  and  arc  denoted  as 
Scheme  1  and  Scheme  2,  respectively. 

The  energy  was  found  to  be  consistent  for  the  two  sets  of  weights  to  <  l.xlO  Harirccs  and  the 
electric  transition  dipole  moments  (for  X'X*  <-  a^Z^  differed  by  less  than  \7c.  An  additional  check 
on  the  choice  of  weighing  factors  comes  from  the  comparison  of  the  computed  molecular  constants 
with  the  experimental  values  (see  Table  4)  when  available.  Finally,  comparison  of  the  results  for  the 
a^Xu  state  from  this  study  with  extensive  non-state  averaged  calculations  of  Konowalow  and 
Lcngsficld  (1987b)  shows  good  agreement  for  the  r^.  D^.,  (o^.,  and  the  description  of  the  "intermediate 
hump"  in  the  potential  energy  curv'c  (PEC)  for  this  slate. 

At  smaller  iniemuclear  separations  (R  =  1.3,  1.5,  1.6),  the  basis  set  bectunc  linearly  dependent  and 
molecular  orbitals  were  eliminated  in  order  to  obtain  convergence  in  the  Cl  diagonalizalion.  At 
R(Hc  -  He)  =  1.30,  one  MO  of  b,„  symmetry  was  eliminated  from  the  virtual  space,  and  at  R  =  1.50 
and  1.60,  two  MOs  of  b,^,  symmetry  were  eliminated.  These  correspond  to  MOs  consisting  primarily 
of  the  most  diffuse  s-type  atomic  orbital  (AOs). 
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The  effect  of  climinaiing  these  MOs  was  checked  at  R  =  1.70  by  comparing  the  results  for 
calculations  with  all  MOs  included  to  calculations  where  first  one  MO  and  then  a  second  MO  was 
removed  (in  decreasing  order  of  diffuseness)  from  the  b,^  IRREP.  It  was  found  that  by  eliminating 
one  MO,  and  then  a  second,  the  energy  differed  by  no  more  ilian  ±3  x  10'^  Hartrce  for  any  state  when 
compared  to  the  calculaiion  using  all  the  MOs.  The  transition  dipole  moment  differed  by  no  more 
than  3%.  The  vibrational  analysis  was  re-run  with  the  electric  transition  moment  increased  at  these 
three  points  by  twice  the  variation  witnessed  at  R  =  1.70  (i.e.,  by  a  factor  of  .06),  then  again  with  the 
transition  moment  decreased  by  a  factor  of  .06  at  these  points.  All  the  resulting  lifetimes  were 
identical  to  the  initial  results  to  within  at  least  two  significant  digits.  Stability  in  the  lifetimes  to  this 
level  of  precision  is  acceptable  for  this  study. 

The  final  zerolh-order  wavefunctions  were  obtained  from  second-order  CIs  with  respect  to  the 
SA-MCSCF  active  space.  The  size  of  the  resulting  Cl  expansion  (in  number  of  CSFs)  for  each  state  is 
X*i;(27,381),  a^i;(38,218),  b^ng(33.702),  F'n,(23,490),  for  the  cases  where  no  MOs  were  eliminated. 
When  one  and  two  MOs  of  bj^  symmetiy'  were  eliminated,  the  corresponding  totals  are  (26,364, 

36,794,  32,607,  22,736)  and  (25,369,  35,403,  31,530,  21.992),  respectively. 

In  the  vibrational  analyses,  PECs  were  represented  by  spline  functions  over  the  region  for  which 
ab  initio  data  was  available  with  extrapolation  using  Lcnnard-Joncs  6-12  functional  forms.  The  a^Z^, 
b^flg,  and  F'n^  states  were  represented  by  spline  functions  for  tltc  region  R  =  1.30  -  6.50  bohr,  wliilc 
spline  functions  were  used  to  represent  the  X'Z^  PEC  for  the  region  R  =  1.30  -  15.0  bohr.  The  total 
electric  dipole  transition  moment  was  also  represented  by  a  spline  function  for  points  along  R  =  1.3  - 
6.5  bohr,  and  described  by  a  second-order  polynomial  outside  this  range. 

4.  RESULTS  AND  DISCUSSION 

4.1  Stale  Propen ies.  One  finds  the  following  state  description  at  R  =  2.00  bohr: 


x‘i;: 

lojio; 

10^20^1 

b^n^: 

loJlCTjjt^ 

F'n  ; 

u 

lo;laj7r^ 
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Much  of  the  behavior  in  the  bound  region  can  be  understood  by  treating  Hej  as  Hej  plus  an 
electron  in  a  Rydberg  orbital.  The  three  electrons  of  He^  form  the  tightly  bound  "core"  electrons  (with 
MO  occupation  lo^lCu)  which  interact  to  form  the  attractive  potential  at  small  R  values.  All  three 
excited  states  have  tliis  core  description  in  the  dominant  CSF  within  the  bound  region.  For  R  >  3.0 
bohr,  the  contribution  from  a  CSF  containing  the  anti-bonding  configuration  lci^,l0u  begins  to  make  a 
significant  contribution  for  the  three  excited  states  in  this  study.  Figure  1  contains  plots  of  the 
potential  energy  curves  (PECs)  for  the  four  stales  of  interest,  and  Table  2  reports  the  actual  energies. 
Table  4  compares  the  spectroscopic  constants  for  the  four  stales  of  interest  as  predicted  by  this  study 
and  experiment.  These  are  provided,  in  part,  as  a  check  on  the  overall  quality  of  the  wavefunctions 
used  in  this  study.  The  llicoretical  values  are  calculated  as  the  difference  in  energy  Efr^)  -  E(R  = 
40  bohr),  with  Efr^)  determined  from  a  three-point  fit  to  a  parabola. 

Table  5  lists  the  lowest  10  vibraiion;il  levels  for  state  as  calculated  from  the  vibrational 
analysis.  The  v  =  9  level  lies  13,332  cm  '  above  the  equilibrium  energy  and  2,848  cm  '  below  the 
barrier  maximum.  Table  5  also  includes  die  v  =  0  -  9  levels  for  the  b^n^  and  the  v  =  0  -  3  for  the 

F'n,. 

In  the  following  sections,  an  analysis  of  the  a^X*,  F'D^,  and  b^FIg  states  of  Hcj  is  presented. 

A  general  discussion  of  the  structure  of  die  wavefunctions  for  the  excited  stales  of  Hcj  and  its 
relationship  to  the  shape  of  the  PECs  can  be  also  found  in  papers  by  MuUikcn  (1964a,  1964b,  1966), 
and  by  Gubennan  and  Goddard  (1975),  who  place  special  emphasis  on  the  X  states. 

4.1.1  The  a^Xu  Slate.  The  small  barrier  to  dissociation,  or  "hump",  has  a  maximum  in  this  study 
at  R  =  2.70A,  and  is  reported  (Jordan,  Siddigui,  and  Siska  1986;  Milliken  1964b)  to  occur  from  the 
competition  between  the  attractive  ionic-like  core  and  the  long-range  repulsive  interaction.  Table  3 
gives  various  estimates  of  this  barrier.  The  present  study  calculates  the  barrier  height  to  be  1.56 
kcal/mol  at  2.70A,  which  agrees  well  with  the  relatively  recent  experimental  value  of  1.43±.05 
kcal/mol  at  2.72±.04A  reported  by  Jordan,  Siddigui,  and  Siska  (1986).  Probably  the  best  llicoretical 
estimate  (and  maybe  the  best  overall  estimate)  for  this  barrier  comes  from  a  recent  paper  by 
Konowalow  and  Lengsfield  (1987)  who  calculate  the  barrier  to  lie  at  2.712A  with  a  height  of  1.507 
kcal/mol.  These  agree  quite  well  with  the  values  obtained  in  the  current  study.  As  can  be  seen  in 
Table  4,  the  calculated  r^,  (o^,  T^,  and  D^.  vary  from  experiment  by  no  more  than  \%. 
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4.1.2  The  p’riy  State.  The  existence  of  the  barrier  with  a  predicted  maximum  of  10.9  kcal/mol  at 
R  =  I.79A  has  been  shown  to  arise  primarily  from  an  avoided  crossing  with  a  higher  state  of  the  same 
^riy  symmetr)',  particularly  the  interaction  with  the  state  that  dissociates  to  He(ls^)  +  He*(ls3d)  (Gupta 
and  Matson  1969).  Mulliken  (1964a,  1966)  aLso  predicted  the  existence  of  a  barrier  due  to  a  change¬ 
over  of  the  p'riy  state  from  one  which  looks  like  a  3dn  state  in  tlie  united  atom  orbital  (UAO) 
description  to  one  with  ls“  •+  ls2p  character  as  it  approaches  Uie  dissociation  limit. 

Table  3  compares  the  current  values  for  the  height  (10.9  kcal/mol)  and  location  (Rmiu^  =  1.79A)  of 
this  barrier  with  results  from  two  other  tlicoretical  studies.  Gupta  and  Matsen  (1969)  calculated  values 
of  13.5  kcal/mol  and  1.73A  for  the  barrier  height  and  location,  while  Browne’s  (1965)  results  predict 
12.5  kcal/mol  at  1.77A.  The  predicted  location  for  the  maximum  from  these  three  studies  arc  in 
reasonable  accord.  It  is  not  surprising  however,  that  the  barrier  height  calculated  in  this  study 
(i.e.,  10.9  kcal/mol)  differs  significantly  from  these  other,  very  early  calculations. 

The  calculated  r^,  co^,  and  T^  for  this  state  are  in  good  accord  with  the  experimental  values  (see 
Table  4).  No  experimental  was  reported.  It  should  be  pointed  out  that  the  b^fl^  and  P'fl^  states  do 
not  enter  independently  into  our  ctilculaiions  of  the  a^Su  transition  moment,  since  the 

contributions  from  a  large  number  of  states  of  a  particular  symmetry  are  obtained  by  solving  for  the 
perturbation  over  CSPs.  However,  it  is  still  important  that  the  CSP  lists  and  MOs  provide  a  suitable 
basis  for  describing  the  ‘fl^j  and  ^Ptg  spaces,  and  so  a  comparison  of  the  tlicoretical  and  the 
experimental  spectroscopic  constants  provides  a  useful  check  of  our  calculations. 

4.1.3  The  b^Hg  State.  The  b^Og  has  a  UAO  description  of  2pn  which  dissociates  to  He(ls“)  + 
Hc*(ls2p),  thus  Mulliken  predicted  that  a  hump  is  not  likely  to  occur  in  the  PEC  for  this  state  since  it 
is  not  a  "promoted"  Rydberg  MO  state.  The  potential  energy  curve  for  the  b^flg  state  shown  in 
Figure  1  docs  not  indicate  a  barrier,  thus  supporting  this  prediction. 

Table  4  shows  the  calculated  r^,  w^.,  and  T^,  to  again  be  in  excellent  agreement  with  the  available 
experimental  values.  No  experimental  is  reported. 
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4.2  Transit  ion  Propcrlics. 


4.2.1  Spin-Orbit  Interactions.  The  first-order  corrections  to  the  x'X*  and  states  arise  from 
interactions  of  these  zerolh-order  wavefunctions  with  the  ^fl^j  and  ’fT^  symmetry  manifolds, 
respectively.  The  magnitude  of  the  perturbation  of  the  by  the  'fly  manifold  is  plotted  in 

Figure  2a,  and  labeled  Curve  A.  Curve  B  in  Figure  2a  represents  the  first-order  SO  perturbation  of  the 
zeroth-order  wavefunction  attributable  to  only  the  lowest  energy  state  of  'fl^  symmetry,  the  F’n^ 
state.  That  is  the  L=1  trunction  of  Equation  5.  Therefore,  the  difference  between  Curves  A  and  B 
reflects  the  error  in  the  first-order  perturbation  treatment  of  4^°(a^Zui)  that  is  being  introduced  by 
truncating  the  summation  in  Equation  5  to  simply  L=l.  The  analogous  information  is  plotted  in 
Figure  2b  for  the  X'Zg  state  being  perturbed  by  the  ^fl^  manifold  (Curve  A)  or  only  the  b^flg  state 
(Curve  B). 

One  can  immediately  sec  that  much  of  the  contribution  to  the  total  perturbation  is  excluded  from 
the  '^’’s  if  only  the  interaction  with  the  lowest  energy  'n^  or  state  is  included.  The  difference  in 
the  contributions  at  R=2.()0  is  a  factor  of  ten  for  the  a^Zuj  -  'fl^j  SO  interaction  and  more  than  a 
factor  of  20  for  tlie  X’Zgo+  -  interaction.  The  discrepancies  change  near  R=4.0  bohr,  where  the 
perturbation  of  the  a’Z*  state  by  a  single  'fl^  state  accounts  for  approximately  79%  of  the  total 
interaction  attributed  to  the  ’0^  manifold.  However,  the  single-state  approximation  for  the  x’Zgo+  - 
perturbation  is  more  than  one  hundred-fold  less  than  that  calculated  from  the  interaction  with  the 
entire  manifold  for  most  of  the  bound  region. 

4.2.2  Electric  Transition  Dipole  Moment  and  Lifetimes.  The  total  electric  transition  dipole 
moment,  l.t,(a^Zu,X*Zp,  obtained  from  the  perturbed  wavefunctions  in  Equations  9,  as  well  as  its 
singlet  and  triplet  components  (as  given  in  Equation  11),  arc  plotted  as  the  dotted  curves  in  Figure  3, 
and  Table  6  lists  the  values  of  l.ii(a^Xu,X'Zp  as  a  function  of  R(Hc-He).  It  can  be  seen  that  the 
singlet  component  dominates  over  most  of  the  a^Zu  bound  potential,  with  the  triplet  component  having 
comparable  magnitude  only  at  small  intcrnuclcar  separations.  At  R=1.6,  the  triplet  component  is 
already  a  factor  of  five  smaller  than  the  single  contribution. 

The  two  moments  have  opposite  signs  for  values  less  than  2.0  bohr,  and  then  have  the  same  sign 
up  through  R=3.5,  where  the  signs  are  once  again  opposites.  The  difference  is  signs  at  small 
intemuclear  separation  cau.ses  a  cancellation  in  forming  the  total  transition  moment,  generating  a  near 
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zero  moment  at  R=1.50  bohr.  From  R>1.85,  the  total  transition  dipole  is  largely  determined  by  the 
singlet  component  which  has  a  maximum  value  of  6.0xl0‘^  au  at  R=3.8  (from  fitting  to  a  parabola). 
The  decreasing  transition  moment  for  large  R  is  consistent  with  the  separated  atom  limit,  for  which  the 
electric  transition  dipole  moment  must  go  to  zero  as  it  represents  a  Hc(^Sg)  — ^  Hc(*Sg)  transition. 

The  single  stale  L=1  approximation  in  Equation  5  is  also  considered  in  Figure  3.  The  solid  curves 
in  Figure  3  provide  tlie  singlet  and  triplet  components,  as  well  as  the  total  PjCa^Su.X'Zp  as  given  in 
Equation  11,  but  calculated  within  the  L=1  approximation.  Comparing  the  dotted  curves  with  the  solid 
curves  one  finds  at  least  three  main  differences.  First,  the  singlet  contribution  to  M.,(a^Su,X^Lp  for  the 
single-stale  perturbation  (SSP)  (solid  curve)  is  essentially  zero  for  the  region  R=1.3  to  2.6  bohr,  in 
sharp  contrast  to  the  singlet  contribution  given  by  the  doited  curve,  which  never  falls  below  50%  of 
the  maximum  in  Pi(a^XpX’Sp.  The  second  observation  is  that  the  triplet  contribution  to 
lti(a^]E;^,X^Sp,  from  the  SSP  is  much  larger  in  this  region.  For  example,  at  1.85  bohr,  the  triplet 
contribution  is  -2.1x10  '’  au  for  the  SSP,  while  it  is  essentially  zero  for  the  perturbation  over  the 
manifold  of  stales.  The  third  feature  is  the  relative  magnitude  of  the  total  |.t,(a^XpX'Lp’s  near  their 
maxima.  For  example,  at  R=4.0  bohr,  we  find  p,(a^XpX‘Zp=6.0xl0  *^  au  for  the  calculation  over  the 
manifold,  while  the  SSP  gives  Pi(a^lpX'Zp=3.5xl0'^  au,  and  therefore  accounts  for  only  58%  of 
the  predicted  total  magnitude  of  the  transition  dipole  moment.  However,  from  Figure  2a,  we  sec  that 
at  this  geometry  approximately  79%  of  the  S-0  perturbation  is  accounted  for  using  the  SSP.  Thus,  the 
electric  transition  dipole  moment  converges  more  slowly  with  respect  to  L,  the  number  of  excited 
states  included  in  Equation  5,  than  the  S-0  first-order  perturbation  contribution  to  'P(a^5^j). 

Table  5  lists  the  predicted  lifetimes  and  energies  from  this  study  for  the  v=0-9  vibrational  levels  of 
the  a^Z*  state  for  a  radiative  decay  proce.ss  to  the  repulsive  X'!*  state.  The  predicted  lifetime  for  the 
v=0  level  is  18  sec,  which  is  consistent  with  the  more  recent  experimental  prediction  of  10  see 
(v=unknown)  for  a  lower  bound  in  liquid  helium.  Tlie  lifetimes  are  seen  to  monotonically  decrease 
with  increasing  vibrational  quantum  number,  at  least  up  to  v=9.  At  v=5,  the  lifetime  falls  below  the 
predicted  lower  bound  of  10  sec.  The  calculated  lifetime  of  the  v=0  level  using  the  electric  transition 
dipole  moment  represented  by  the  solid  curve  in  Figure  3  (from  the  single  stale  approximation  to 
Equation  5),  is  predicted  to  be  195  sec,  in  sharp  contrast  with  results  determined  by  including  all  of 
the  eigenstates  in  our  CSF  basis. 
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5.  CONCLUSIONS 


The  lifetime  for  the  Hcj  excited  state  is  predicted  to  be  18  sec  for  the  v=0  vibrational  level  in 
the  gas  phase,  supporting  the  experimental  value  for  the  lower  bound  (in  condensed  phase)  offered  by 
Mehrotra,  ct  al.  (1979),  of  10  sec.  These  calculations  also  predict  the  lifetime  to  decrease  continuously 
with  increasing  vibrational  quantum  number,  at  least  up  to  tlte  v=9  vibrational  state. 

One  finds  that  the  Pi(a^Xu.X'Xp  shows  maxima  near  4  bohr,  and  the  electric  transition  dipole 
moment  for  intemuclear  separations  greater  titan  1.60  bohr  is  determined  almost  entirely  by  the  singlet 
component,  <4''(’nui:a^X*i)llx^,l'F°(X'X^yJ>.  S-O  interactions  originating  in  'fl^  states  beyond  the 
F'fly  arc  essential  to  the  characterization  of  the  'F'('ny,;a^X„i)  wavefunction,  as  well  as  p.,(a^Xu,X^Xp. 
This  is  a  strong  argument  in  favor  of  using  the  method  employed  in  this  study,  which  is  designed 
specifically  to  include  these  higher  energy  contributions  at  little  or  no  additional  cost. 
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Figure  lb.  Blow-up  oF  the  Potential  Energy  Curves  for  the  a^Z...  b^H^. 


and  F*n..  States. 


S-O  Perturbation  xl0(+6) 


R(He-He)  lau) 


F'gu''c2a.  First-Order  Contribution  to  the  Spin-Orbit  Pcnurhation  of  the  hv  the  *n..  State 
Manifold  CCurve  A)  and  by  the  F'n,.  State  tCurve  B). 
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S-O  Perturbation  x  10 1+61 


RlHe-Hel  !au) 


Figure  2b.  First-Order  Contribution  to  the  Spin-Orbit  Perturbation  of  the  bv  the  ^FI^,  State 
Manifold  tCurve  A)  and  by  the  b^Ff^  State  (Curve  B). 
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Figure  3.  Total  Electric  Dipole  Transition  Moment  and  Singlet  and  Triplet  Components  Using 
Perturbed  Wavefunctions  Including  the  SO  Interactions  with  the  Entire  'O..  and  State 
Manifolds,  f— V,  Electric  Pinole  Transition  Moment  and  ComtxinenLs  for  the  Single-state 
Perturbed  Wavefunctions  Havine  SO  Interactions  With  Only  the  F‘n^,.  and  b^n^- 


17 


Table  1.  Atomic  Basis  Set 


Type 

Exponent 

Coefficient 

s 

501.5045 

0.002498 

75.31147 

0.019099 

17.20769 

0.092978 

4.886925 

0.311074 

1.569584 

1.0 

0.541551 

1.0 

0.193932 

1.0 

0.104560 

1.0 

0.026725 

1.0 

0.008017 

1.0 

P 

10.19643 

0.092050 

2.414857 

0.474058 

0.746691 

1.0 

0.139276 

1.0 

0.032392 

1.0 

0.012 

1.0 

d 

1.5 

1.0 

0.042 

1.0 
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Table  2.  Cl  State  Energies  as  a  Function  of  R(He-He)“ 


r(He-He) 

F‘n„ 

1.30 

-5.302222 

-5.021207 

-4.986410 

-4.901985 

1.50 

-5.469882 

-5.098117 

-5.067915 

-4.986402 

1.60 

-5.532127 

-5.119922 

-5.091493 

-5.011180 

1.70 

-5.583158 

-5.134085 

-5.107148 

-5.027931 

1.85 

-5.642728 

-5.144880 

-5.119785 

-5.042067 

1.90 

-5.658859 

-5.146409 

-5.121831 

-5.044584 

2.00 

-5.686530 

-5.147189 

-5.123516 

-5.047177 

2.10 

-5.709037 

-5.145601 

-5.122740 

-5.047252 

2.15 

-5.718626 

-5.144177 

-5.121637 

-5.046572 

2.30 

-5.741993 

-5.138132 

-5.116346 

-5.042518 

2.40 

-5.753907 

-5.133160 

-5.111710 

-5.038691 

2.50 

-5.763483 

-5.127971 

-5.106540 

-5.034360 

2.60 

-5.771267 

-5.122332 

-5.101183 

-5.029804 

2.70 

-5.177500 

-5.166891 

-5.095719 

-5.025212 

2.85 

-5.784629 

-5.109085 

-5.087677 

-5.018643 

3.00 

-5.789709 

-5.102519 

-5.080019 

-5.012907 

3.25 

-5.795223 

-5.092835 

-5.068858 

-5.006811 

3.40 

— 

— 

-5.005982 

3.50 

-5.798365 

-5.085661 

-5.059869 

-5.006494 

3.70 

-5.799860 

-5.081535 

-5.054116 

-5.008882 

4.00 

-5.801126 

-5.077456 

-5.047474 

-5.013024 

4.25 

-5.801675 

-5.075460 

-5.043445 

-5.015822 

4.50 

-5.801973 

-5.074317 

-5.040475 

-5.017919 

4.75 

-5.802132 

-5.073739 

-5.038312 

-5.019440 

5.00 

-5.802214 

-5.073526 

-5.036750 

-5.020526 

5.25 

-5.802254 

-5.073542 

-5.035628 

-5.021295 

5.50 

-5.802272 

-5.073692 

-5.034824 

-5.021838 

6.00 

-5.802279 

-5.074169 

-5.033841 

-5.022489 

6.50 

-5.802275 

-5.074671 

-5.033337 

-5.022815 

10.00 

-5.802260 

-5.075958 

-5.032746 

-5.023218 

40.00 

-5.802255 

-5.075983 

-5.032699 

-5.023338 

*  Atomic  units  used  throughout. 
Calculated  only  the  F^Fl^. 
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Table  3.  Barrier  Heights  and  Barrier  Positions  for  the  and  F*^^  States 


*  Energies  in  kcal/mol. 

*’  Distances  in  Angstroms. 

'  Peach  (1978). 

Sunil  et  al.  (1983),  MCSCF  calculations. 

‘  Lundlum,  Larson,  and  Caffrey  (1967) 

^  Brutschy  and  Haberland  (1979). 

*  Jordan,  Siddiqui,  and  Siska  (1986). 

Gupta  and  Matsen  (1969),  Valence-bond  calculations. 

‘  Browne  (1965)  did  not  report  a  barrier  position  from  any  fitting  procedure,  so  we  calculated  the  position  and  height 
by  fitting  the  pmteniial  energy  data  in  Table  1  of  Browne  (1965)  to  a  parabola  giving  these  results. 

^  Large-scale  MCSCF  plus  second-order  Cl  (Konowalow  and  Lengsfield  1987). 
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Table  4.  Molecular  Constants  for  the  a^5^,  b^rig,  and  F‘1\  Electronic  States* 


Property 

a'l: 

b^n^ 

F% 

r^  Theory 

1.0493 

1.0681 

1.0869 

Exp. 

1.0457 

1.0635 

1.0849 

rp  b 

^  e 

143,768. 

148,962. 

165,665. 

144,048. 

148,835. 

165,971. 

(5,194.) 

(21,897.) 

(4,787.) 

(21,923.) 

< 

1,816. 

1,766.  . 

1,673. 

1,809. 

1,769. 

1,671. 

15,636. 

19,942. 

5,293. 

15,806. 

*  All  distances  in  angstroms  and  energies  in  cm-1.  Experimental  data  from  Huber  and  Herzberg  (1979). 

The  first  set  of  values  are  with  respect  to  the  at  R=40  au,  and  the  parenthetical  values  are  T^’s  with  respect 
to  the  Ej  of  a’EJ|, 

'  Theoretical  (ojs  from  AG(2-1)  -  AG(l-0)=-2cOjXj  and  (£)^=G(l-0)  +  See  Herzberg  (1950),  pg.  95. 

‘‘  Determined  from  the  energy  difference  between  r^  and  R=:40  au. 


Table  5.  Results  from  Vibrational  Analyses  of  the  b^Og,  and  F*^^  States  with  Energies  in  cm  ' 
and  Lifetimes,  z,  in  Seconds 


V 

Energy 

T 

b^Og 

Energy 

F'n 

Energy 

0 

899 

18 

873 

826 

1 

2,635 

15 

2,570 

2,420 

2 

4,290 

13 

4,199 

3,936 

3 

5,867 

12 

5,757 

5,270 

4 

7,373 

11 

7,242 

5 

8,785 

9.6 

8,658 

6 

10,097 

8.5 

10,000 

7 

11,306 

7.5 

11,270 

8 

12,433 

6.7 

12,459 

9 

13,452 

6.0 

13,569 
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Table  6.  The  Total  Electric  Transition  Dipole  Moment  iai(a^I^,X*Zp,  for  a^SiJ  — »  X'Zg  as  a  Function 
R(He-He)  (in  atomic  units). 


R(He-He) 

<x'z;i^„la'i;> 

(xl0**6)  au 

1.30 

-3.276 

1.50 

1.030 

1.60 

2.046 

1.70 

2.906 

1.85 

3.622 

2.00 

4.058 

2.10 

4,269 

2.15 

4.352 

2.30 

4.557 

2.40 

4.670 

2.50 

4.794 

2.60 

4.868 

2.70 

4.959 

2.85 

5.089 

3.00 

5.340 

3.25 

5.615 

3.50 

5.877 

4.00 

5.990 

4.25 

5.768 

4.50 

5.383 

4.75 

4.882 

5.00 

4.319 

5.50 

3.182 

6.00 

2.199 

6.50 

1.466 

10.00 

2.145(-3)‘ 

40.00 

5.400(-5) 

*  Characteristic  base  ten  noted  parenthetically. 
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IM  ENTIONALLY  LEFT  HLANK. 
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2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 
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